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Abstract 

Quantum fluctuations of scalar fields during inflation could determine the very 
large-scale structure of the universe. In the case of general scalar-tensor gravity 
theories these fluctuations lead to the diffusion of fundamental constants like the 
Planck mass and the effective Brans-Dicke parameter, lo. In the particular case 
of Brans-Dicke gravity, where lo is constant, this leads to runaway solutions with 
infinitely large values of the Planck mass. However, in a theory with variable to 
we find stationary probability distributions with a finite value of the Planck mass 
peaked at exponentially large values of lo after inflation. We conclude that general 
relativity is an attractor during the quantum diffusion of the fields. 
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I. INTRODUCTION 



One of the most important problems for cosmology is the issue of initial conditions for the 
Big Bang The first models of inflation assumed that the universe started in a very hot state 
that supercooled in a metastable vacuum, which then decayed to the true vacuum through 
a first-order phase transition or just rolled down through a second-order phase transition g. 
Chaotic inflation Q opened up the possibility of starting inflation from a wide range of initial 
conditions, including the Planck scale. 

Quantum fluctuations in the inflaton field could produce the small perturbations observed in 
the background radiation §J|, from which galaxies later evolved. Once we include quantum 
fluctuations of the scalar fields during inflation, we find that these can be large and dominate 
the classical evolution as we approach the Planck scale Q. As a consequence, the scalar fields 
diffuse, very much like a particle in Brownian motion. The universe is then divided into causally 
independent inflationary domains, in which the fields acquire different values. One of the most 
fascinating features of inflation is the process of self-reproduction of the universe JO , by which 
the values of the fields in some inflationary domains diffuse towards larger rates of expansion, 
producing new domains and so on for ever. This process might still be occurring, at scales much 
larger than our present horizon. The global behavior of the universe can then be described with 
the formalism of stochastic inflation, using probability distributions for the values of the fields 
in physical space. 

It has recently been shown that there are stationary solutions for the diffusion of the inflaton 
in general relativity A natural question to ask is whether this picture is still valid as 
we approach the Planck era, where quantum fluctuations of the metric become important 
H. Although it is generally assumed that the dynamics of the universe can be described by 
general relativity, the effective theory of gravity might be very different close to the Planck 
scale. So far the only consistent but by no means definite, since we lack the experimental 
observations needed to confirm it, theory of quantum gravity is string theory [||. String theory 
contains in its massless gravitational sector a dilaton scalar field as well as the graviton. The 
low-energy effective theory from strings has the form a scalar-tensor theory of gravity |7j] , with 
non-trivial couplings of the dilaton to matter Therefore, it is expected that the description 
of gravitational phenomena, and in particular inflation, close to the Planck scale should also 
contain this scalar field [ |l(^fif| . 

The string dilaton field can be understood as a Brans-Dicke field Jl2| |, which acts like a 
dynamical gravitational 'constant'. Jordan-Brans-Dicke theory is the simplest scalar-tensor 
theory, with a constant kinetic coupling u>, which is bounded by primordial nucleosynthesis 
|f3| and post-Newtonian experiments to be u> > 500. String theory predicts u> = — 1 in ten 
dimensions ||. However, the low energy effective value of u> depends on the unknown details of 
the compactification mechanism and supersymmetry breaking lq] . In general, one would expect 
a functional dependence of u> on the dilaton field || . Such models were proposed in jl4],[l5| for 
solving the graceful exit problem of extended inflation, and later suggested to be the generic 
asymptotic behavior of scalar-tensor theories' approach to general relativity during the matter 
dominated era ]l6| , |T7| | . 

In this paper we study the very large scale structure of the universe assuming that the 
gravitational interaction is described during inflation, close to Planck scale, by a general scalar- 
tensor theory of gravity, in the context of the stochastic inflation formalism. In Brans-Dicke 
stochastic inflation jT^|il|, we found runaway solutions to the diffusion of the dilaton and 
inflaton fields due to the fact that the Planck boundary, for generic chaotic potentials, is a 
line and the probability distribution does not become stationary, as occurs in general relativity 
||], but slides along this boundary. As a consequence, the value of Planck mass at the end of 
inflation is not well defined (it would depend on new dynamics at large values of the fields, 
e.g. quantum loop corrections to the potential pc|). We will study a particular case in which 
the Brans-Dicke parameter has a simple pole with respect to the Brans-Dicke field, which will 
actually give stationary probability distributions for the diffusion of the inflaton and dilaton 
fields in physical space. We find a probability distribution peaked about domains that produce 
general relativity as the effective theory of gravity at late times and therefore we conclude that 
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it is most probable to live one of those domains .0 This prediction seems to be in good agreement 
with observations. 



II. RUNAWAY SOLUTIONS IN BRANS— DICKE COSMOLOGY 

In this section we review the problem of runaway solutions in Brans-Dicke models of stochas- 
tic inflation. For a detailed analysis see Refs. |l8|,[l9|]. Let us consider the evolution of the inflaton 
field a with a generic chaotic potential in a JBD theory of gravity with dilaton field 4>, 



S = 



d A x, 



1 



1 
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(2.1) 



where Planck mass is written in terms of the dilaton field as Mp(ip) = ^j(f> 2 - For generic inflaton 
potentials of the type V{a) — Xa 2n /2n, the equations of motion for the homogeneous fields in 
the slow-roll approximation are |22j 
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They correspond to circular motion along the lines <fi 2 + ^a 2 — constant, as shown in Fig. 1, 
with the angular variable z = \Jn/2(<fi/a) evolving as Jl9| 



HZ , n\ 

— I + z 2 

UJ 



(2.3) 



The end of inflation occurs for \H\ = H 2 , corresponding to the line z 2 ~ w/2, where the 
slow-roll approximation for the inflaton field breaks down, while the slow-roll approximation 
for the dilaton simply requires u> ^> 1. The Planck boundary is the curve V(a) ~ Mp ((/>), 
4> 2 ~ (cj\/A/27rv / 2n) a n . In the simplest case n = 2, it is the line z 2 = luVX/Ait. 

Apart from classical motion there are quantum fluctuations that act on the background fields 
as stochastic forces and produce a random motion of those fields, very much like Brownian 
motion. The amplitude of quantum fluctuations of S(f> and 8a, whose wavelengths are stretched 
beyond the horizon, can be computed in the slow-roll limit, as in [jl9 23 1, by solving the linearized 
perturbation equations in a de Sitter background (R = 12H 2 = const). They turn out to be 



50 = (4TTk 3 \u k \ 2 ) 
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5a = (4nk 3 \v k \ 2 ) 1/2 = g 
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(2.4) 



Quantum fluctuations of the fields act as stochastic forces on the classical background fields. 
We will describe the quantum diffusion of the coarse-grained fields in terms of the probability 
distribution P p (a 1 <fi; t) to find, at a given time t in a domain with a given physical volume, the 
fields with mean values 6 and a fl8|, 
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(2.5) 



3HP V 



*Note that Coleman's mechanism for the vanishing of the cosmological constant in the context of 
Brans-Dicke theory also predicts general relativity as a low energy effective theory of gravity |21J]. 
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where the first and second terms in each bracket correspond to the classical drift and quantum 
diffusion, respectively. The last term takes into account the different quasi-exponential growth 
of the proper physical volume in different parts of the universe. Those few domains that jump 
in the opposite direction to the classical trajectory contribute with a larger physical space 
and therefore dominate the proper physical volume of the universe. Such domains will split 
into smaller domains, some with lower values of the scalar fields, some with higher. As a 
consequence of the diffusion process, there will always be domains which are still inflating, and 
this corresponds to what is known as the self-reproduction of the inflationary universe jjj . 

Beyond a certain point, the quantum fluctuations Sz of the fields dominate their classical 
motion in a time interval At = H^ 1 and the universe enters the self-reproduction regime. In 
Brans-Dicke inflation this corresponds to the line z = z s , such that 

Sz = zH- 1 . (2.6) 

For n — 2, this is a straight line with ~ oj 3 X/12tt 2 <C 1. For z s < z < z p , new domains create 
more inflationary domains and so on for ever. (The requirement that z s < z p is guaranteed 
by the slow-roll condition, lu 3> 1.) In this picture, most of the volume of the universe today 
is occupied by regions that are still inflating close to Planck boundary, while our own causal 
domain is thought to have evolved from one of those inflationary domains, down the inflaton 
potential, through reheating and into the radiation and matter dominated eras. 

The probability of finding given values of the scalar fields in a given domain in physical 
space can be computed, in a first approximation, with ordinary diffusion equations. In general 
relativity, where the Planck boundary corresponds to a certain value of the inflaton field, it 
is possible to find stationary solutions to the quantum diffusion of the inflaton Q]. However, 
in Brans-Dicke gravity, diffusion occurs in the two-dimensional space (ct, 4>), where the Planck 
boundary is a line. For the usual chaotic inflation potentials, the diffusion of the probability 
distribution along the Planck boundary is unbounded and the fields evolve indefinitely along 
this line [jl8|,[l9| towards larger and larger values producing what we called runaway solutions. 
The diffusion equation for the probability distribution P p (<p,a,t) along the Planck boundary, 
where Hp = (8n/3)Mp((f>) = (8Tr/3)V(a) 1 ^ 2 and diffusion dominates classical drift, can be 
written as 

There are no stationary solutions for generic chaotic potentials of the type V(a) ~ Act 2 ™. It 
was argued in Ref. jl8| that a possible way to stop this runaway behavior of the fields is to 
introduce a cutoff in the potential at ct = <7b due, for instance, to quantum corrections. It is 
also possible to find stationary solutions within Brans-Dicke stochastic inflation by considering 
a non minimal coupling of the inflaton to the curvature scalar, see Ref. |2(i| ] . 

In this paper we explore another possibility which arises in more general scalar-tensor grav- 
ity theories where the Brans-Dicke parameter lu becomes a function of the dilaton field. Such 
theories possess the observationally desirable feature of including the general relativistic be- 
havior in the limit lu — ► oo. As we shall see, these models can also yield stationary probability 
distributions. 



III. VARIABLE OMEGA PARAMETER 



Let us now consider the classical evolution of the inflaton field with a generic chaotic potential, 
in the context of an arbitrary scalar-tensor theory of gravity, 
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Here the parameter w, a constant in Brans-Dicke theory, becomes a function of the effective 
Planck mass 
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M 2 (0) = 16tt/(0) , 

m (3 - 2) 



2[/W 



where 167r/(0) acts like the Brans-Dicke scalar JE^ |. The equations of motion of theory ( |3.l| ) 
can be written as 

v 2 = f{4>)R, 

(3.3) 

V 2 a - - V'{a) , 



2/(0) ( Rpu - \g^R ) = 9»vV{a) + 2 (V M V„ - 9lut V 2 ) /(0) 



1 A / 1 ( 3 - 4 ) 

d^d v (j) - -g^(d(j)) 2 J + f d^adya - -g^{da) 2 



We can then write the exact equations for the homogeneous fields in a spatially flat (k = 0) 
Friedmann-Robertson- Walker metric as 

(2w + 3) (/(0) + 3#/(0)) + co'(f)f^) 2 = 2V(a) - \a 2 , 

v ' (3-5) 

a + 3H& = - V'(a) , 



/ ( - .) ( 6H + 12H 2 ) = 2V(a) - 3/(0) - 9#/(0) - X -a 2 - \^ 2 , 

1 l . (3-6) 
/(0) QH 2 = V(a) - 6iT/(0) + ^ + ^ 2 , 



where H = a/a, and a is the scale factor. Note that for /((/>) = <t> 2 /8u> we recover the usual BD 
equations. 

During inflation, we can write the equations of motion of the homogeneous fields and <r, in 
the slow-roll approximation, as 



1/((T) 47T OCT 



ff2 = W 



6/(0) ' 

For the slow-roll solution to be an attractor, /(0) must satisfy the following conditions, 

/'M 2 « /«.) =* A!) » i 

v (f) 

In addition we have the straightforward generalization to scalar-tensor gravity of the familiar 
slow-roll conditions for the inflaton potential, f{4>){V' /V) 2 < 1 and f{<j))\V"/V\ < 1. 

The end of inflation in our theory occurs when \H\ = H 2 , or V — a 2 + <f> 2 + 3/(0) + 3Hf(4>). 
For the generic chaotic potential V(a) — \ Ac 4 , the end of inflation corresponds, in the slow- 
roll approximation, to a e ~ Mp(0 e )/ yft. Note that in our theory, inflation ends when the 
inflaton field starts oscillating around the minimum of its potential, while the dilaton field 
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becomes essentially constant. On the other hand, the Planck boundary is approximately given 
by V(a) = M${4>) §, or a p ~ V% \~ 1/4 M P (q> p ). 

The quantum fluctuations in a variable oj(f) model remain those of Eqs. ( |2.4| ) since they 
were calculated in the slow-roll limit, which should still hold here. The variable u>(f) then only 
affects how the diffusion of <f> is reflected in f(4>), and thus in the Hubble rate H. Along the 
Planck boundary the rate of expansion is given by Hp = (87r/3)ATp(0) and thus will still diffuse 
towards larger values of Mp. 

A simple form to consider is uj(f) = ujq + uj m f m 0], which satisfies the slow-roll conditions 
for o^o 3> 1 and u) m > 0. For m > we expect that quantum fluctuations of the fields will 
drive the effective value of the Brans-Dicke parameter to infinity, recovering general relativity 
in this limit. However, in general, we will not recover a finite value of Mp, i.e. there will still 
be runaway solutions for the quantum diffusion. A divergent lu is not sufficient to recover a 
stationary probability distribution. However, in the next section we analyze a scalar-tensor 
theory that has an upper bound on the Planck mass which does produce both a stationary 
distribution and a divergent u). 



IV. THE MODEL 



In this section we study the classical evolution of the theory fl3.1| ) , for a particular non minimal 
coupling f(<p). We are interested in functions that acquire a maximum (/'(</>) — 0) at some 
value of 4>, which will lead to a singular value of Lu(f), see Eq. (|3.2|). We choose a simple case 
for which there are analytical solutions, 



f{4>) = fo sin 2 acj) . 
In this case, the Brans-Dicke parameter behaves as, 

UJ LOO 



"(/) 



1 - f/fo 



(4.1) 



(4.2) 



where a 2 = (8wo/o) _1 = ZkGq/ojq. Equation (4.2) has a simple pole at / = /o, where we 
recover general relativity with agravitational constant Go- This type of pole behavior for the 
lu parameter was used in Ref. |l5| to allow for a graceful exit of extended inflation. Note that 
it behaves as ordinary Brans-Dicke t heor y with lu = u>q for a(f> <C 1. 

The slow-roll equations of motion (3.7) for the theory Via) = A<7 4 /4 become, 
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3i0r 
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H = 
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V 3 / sin a 
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(4.3) 
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We can redefine u 
become 



{X/Sloq) 1 / 2 t/a , x = ao , y = sina</> , under which the equations of motion 



y' = x 2 (l~y 2 ), 
x' = -xy , 



(4.4) 



where a pri me here denotes a derivative with respect to u. The slow-roll conditions for the 

■ ;1 and luq^> y 2 , which are both satisfied by ujq ^ 1- 

• 1 '"' 1 y 2 ) — constant, see 



dilaton (3 



require ujq 1 — y" 
There are exact solutions to the system of equations (^4), x 2 — ln(l 
Fig . 1. For small values of y, we recover the circular solutions found in the Brans-Dicke case 
|l9[ , while in general we have 



G 



= exp (a;f - x 2 ) 



(4.5) 



Note that during inflation ui increases exponentially. This behavior will be important for the 
approach to general relativity. 

In the variable z = y/x, the end of inflation and the Planck boundary correspond to[j] 

— 2 ' 

/I ( 46 ) 

2 ^0V A 
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The simple analytic results for slow-roll inflation in this model allows us to describe very simply 
the evolution of various quantities such as the number of e-foldings from the end of inflation, 



N = ujq In 



y e ( l -y 



v 2 V 1 - vl 



(4.7) 



which approaches N ~ ujq x 2 for x 3> x e . 

The new feature in this particular model is the existence of a boundary in the motion of the 
dilaton. As it grows towards y = 1 it slows down and eventually stops while the inflaton field 
x ends inflation. This effect will be critical for the behavior of the probability distribution, as 
we will see shortly. It is simply a consequence of the local maximum in the function /(</>), and 
will occur in a wide range of functions. 



V. SELF-REPRODUCTION AND STATIONARY DISTRIBUTIONS 

We will now study the onset of the self-reproduction of the inflationary universe in our model. 
The bifurcation line is defined by the values of the scalar fields for which the maximum of the 
probability distribution P p starts increasing, or equivalently, where the quantum diffusion of 
the fields dominates its classical motion in the time interval H . The bifurcation line can be 
written, in the z variable, as 



There will be self-reproduction of the universe for z s > z > z p . The probability distribution 



3 \ 1/2 



P p (4>, &,t) follows the diffusion equation (2.5), although, as in the Brans-Dicke case discussed in 
section jj, it will very rapidly diffuse along the angular direction z, towards the Planck boundary, 
which is an absorbing boundary It is possible to estimate the relative dispersion in z 

due to quantum fluctuations, 

7 = A^- y2) + z ^ /2 ■ (5 - 2) 

Z ^QZ 

Note that the prefactor z 2 /uiq z 2 < 2/y/3ujQ for z > z p so the dispersion must be small for small 



values of z and decrease as we approach y = 1, which strongly suggests that we can approximate 
the diffusion as being essentially one dimensional, along the Planck boundary. The distribution 
then evolves along this boundary following Eq. ( |2.7P , towards large values of a until it reaches 
the line y = 1, which acts like an effective cutoff and the distribution acquires a peak. Diff usio n 



in (f> across the boundary y = 1 is automatically identified with values of y < 1, see Eq. flf.l|) 



The upper limit on a is then a = (2/GovA) 1 / 2 = at,. At this point, the parameter cj(/) becomes 



^Note that these are the same values as in Brans-Dicke theory 
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infinite, while foj'(f)/w 2 — » 0, corresponding to the general relativistic post-Newtonian limit 



Let us compute the shape of the distribution close to the peak, since that will give us some idea 
of possible deviations from general relativity. The diffusion e quat ion ( |2.5| ) for the probability 
distribution in the physical frame, along the Planck boundary (2.7), for z p <C 1, can be written 
in the form of a Schrodinger equation, 



d 2 ^ 
P p {<j,t) 



W(s)^(s 
= exp 



(5.3) 



where 'J'(s) is the wave function, s 2 = 1/cr, and W(s) = — 72n/s 2 ^/\ is the singular potential. 
It is well known from quantum mechanics that there are no regular solutions of (5.3) at s = 
for a singular potential like — 1/s 2 . However, in our case, the boundary at a — <7& permits the 
existence of stationary solutions, with ^(crf,) = 0. A numerical solution is presented in Fig. 2. 
Unfortunately, we do not have an analytic expression for it, but one can easily write down 
WKB approximate solutions as 



*(cr) ~ ( 1 | exp 




a > er* , 



(5.4) 

where a* = eV\/72tt = a b (l - (9tt\/A/128) 1/3 ) is the WKB turning point. This solution has 
a very sharp maximum close to ct* and an exponential decay for small a. 

It is now possible to study the way general relativity is approached as the probability distri- 
bution moves towards the critical point <Jb- The value of the variable Brans-Dicke parameter 
along the Planck boundary can be written as 



cj 



1 



(5.5) 



which shows a pole singularity at <Jb- Thus near the peak of the probability distribution function 
we find 



wp(cr*) 



1/3 



Wo 



(5.6) 



which may be quite large for the values of A ~ 10~ 12 required by density perturbations 
However, the main increase of to comes from the classical evolution of those domains that 
started close to the peak of the distribution and later evolved towards the end of inflation.^ It 
is possible to compute this increase using the slow-roll approximate solution (4.5). By the end 
of inflation, the effective value of ui coming from domains at the peak of the distribution takes 
the expression 



7 end 



wp(cr*) exp 



2 /2tt ^ 



(5.7) 



"Note that diffusion is important between the Planck boundary and the self-reproduction boundary, 
however one expects that, as in the case of Brans-Dicke j^], the peak of the distribution will closely 
follow the classical trajectory. 
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which is exponentially large, compared to c^o- We conclude that stochastic inflation drives the 
effective value of u> at the end of inflation to be exponentially large, with a sharply peaked 
probability distribution. It is important to emphasize here that once the distribution on the 
Planck boundary is peaked at er» , it is the classical motion towards the end of inflation which 
is responsible for the exponential increase of the effective Brans-Dicke parameter. 

VI. CONCLUSIONS 

Why do the constants of nature take the values we observe? Couplings range fourteen orders 
of magnitude; masses are smaller than 10 -17 Afp and range eleven orders of magnitude; the 
vacuum energy is smaller than 10 _120 Mp, and so on. A possible answer, very popular among 
particle physicists, is that there is a unique logically consistent theory of everything, where all 
fundamental constants are determined from its vacuum state. Unfortunately, this state is prob- 
ably not unique, e.g. in superstrings it strongly depends on the compactification mechanism 
and supersymmetry breaking 0]. On the other hand, quantum cosmology proposes that the 
so-called wave function of the universe provides a probability distribution for all fundamental 
constants. It is usually studied in the canonical or Euclidean approach, which has problems of 
interpretation related to the choice of measure. In Ref. § it was suggested that the stochastic 
inflation formalism could provide a reasonable framework within which to answer these ques- 
tions. Here an exponentially large, causally disconnected inflationary domain replaces a single 
nucleated universe of Euclidean quantum cosmology. This formalism proposes that the global 
measure should be given by the probability distribution in physical space [^[^H), which takes 
into account the proper volume of the universe. 

Stochastic inflation describes the quantum diffusion of fields close to the Planck boundary. 
It uses branching diffusion equations to derive the probability of finding a given value of the 
scalar fields that drive inflation in a given physical proper volume [Q. It can be analyzed 
in the context of general relativity or in other theories of gravity, like scalar-tensor theories, 
where the gravitational coupling (Newton's constant in general relativity), becomes another 
dynamical field, the Brans-Dicke field. In the stochastic picture this leads to different values 
of the effective Planck mass in different exponentially large, causally disconnected, parts of 
the universe. This picture may be incorporated in a theory of evolution of the universe ^6j, 
where quantum fluctuations of Planck mass could act as a mechanism for mutation, while a 
selection mechanism establishes that its value should be as large as possible in order to increase 
the rate of expansion, and therefore the proper volume of the universe. Unfortunately, in 
the simplest scalar-tensor theory, Brans-Dicke theory, diffusion close to the Planck boundary 
leads to runaway solutions where the global volume of the universe becomes dominated by 
regions with infinitely large Planck mass, in conflict with observations unless new dynamics is 
introduced into the model Jl8| p0| . 

In scalar-tensor theories with power-law behavior of the Brans-Dicke parameter, there are 
still runaway solutions, but by considering a scalar-tensor theory with an upper bound on the 
BD field (corresponding to a pole in the BD parameter), we have shown that not only do we 
recover a stationary probability distribution for the fields along the Planck boundary (peaked at 
the maximum allowed value of the Planck mass) , but that the low-energy effective cu parameter 
becomes exponentially large, thus recovering the general relativistic behavior. It is important 
to emphasize that this result is expected to be generic in all models involving a maximum value 
of the Planck mass, not just the particular model considered here. 

This purely quantum diffusion process towards large values of to is then reinforced by the sub- 
sequent classical evolution of the inflationary universe during which the Brans-Dicke parameter 
exponentially approaches the general relativistic limit. The ability of infinite to to act as an 
attractor in classical cosmology is well-known during inflation fl5|| and matter dominated era 
|l6| . What we have presented in this paper is a quantum process which enables us to attribute 
a relative probability of finding a given value of lo in the post-inflationary universe. Along with 
the classical evolution, this quantum diffusion mechanism predicts an effective theory of gravity 
which at latc-timcs is indistinguishable from general relativity. 
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FIGURE CAPTIONS 



FIG. 1. The classical evolution of dilaton and inflaton fields during inflation in the (x,y) plane, 
see Eq. (4.4), for the scalar-tensor theory defined by (4.1) is represented by the continuous curves, 
while the dashed curves correspond to the solutions in ordinary Brans-Dicke theory. Classical motion 
starts at Planck boundary (z p ) and ends at the end of inflation boundary (z e ), represented by the thick 
straight lines, while the dotted line corresponds to the self- reproduction boundary (z s ). The dot-dashed 
horizontal line corresponds to the general relativistic limit y — » 1. Note that the classical motion in 
our theory very quickly approaches that limit. 

FIG. 2. The probability distribution 9(a) along the Planck boundary, for at = 1 and a, = 0.812. 
The WKB solutions (5.4) are good approximations to the numerical result away from the turning point 
<r» , marked here with a vertical line. 
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Fig.2 




